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, Abstract. We have observed the flare star YYGem simultaneously 

with XMM-Newton and Chandra as part of a multi-wavelength campaign 
t^j- , aiming at a study of variability related to magnetic activity in this short- 

CK ' period eclipsing binary. Here we report on the first results from the 

analysis of the X-ray spectrum. The vicinity of the star provides high 
enough S/N in the CCD cameras onboard XMM-Newton to allow for 
time-resolved spectroscopy. Since the data are acquired simultaneously 
they allow for a cross-calibration check of the performance of the XMM- 
Newton RGS and the LETGS on Chandra. 
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1. Introduction 

x 

■ YY Gem is the optically faintest of the three visual binaries in the Castor sextu- 

plet. It is itself an eclipsing spectroscopic binary with period of 0.81 d. The two 
components of YY Gem are both of spectral type dMle, and belong to the class 
of BYDra variables. Indeed, YYGem was the first stellar system on which pe- 
riodic photometric variability was detected (Kron 1952). Since the discovery of 
X-ray emission from the Castor system by the Einstein satellite, the system was 
studied by virtually all X-ray observatories (Vaiana et al. 1981, Pallavicini et al. 
1990, Gotthelf et al. 1994, Schmitt et al. 1994, Giidel et al. 2001). Flares on 
YY Gem have been recorded from all parts of the electromagnetic spectrum. The 
extraordinary activity of this object may be related to its binarity (the frequency 
of photometric flares seems to be enhanced in the interbinary space suggesting 
interaction between the magnetospheres of the two stellar components; Doyle Ik 
Mathioudakis 1990), and makes it a prime target for simultaneous monitoring 
at different wavelengths. 
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2. Observations and Data Analysis 

YYGem was observed by both Chandra and XMM-Newton on Sep 29/30, 2000 
for a total observing time of 59 ksec and 55 ksec, respectively. The XMM-Newton 
observations were obtained in the full-frame mode of EPIC-pn, with the thick 
filter inserted for both pn and MOS. We perform the data analysis with the 
standard XMM-Newton Science Analysis System (SAS). Chandra was used in 
the LETGS configuration, i.e. the Low Energy Transmission Grating (LETG) 
combined with the High Resolution Camera for Spectroscopy (HRC-S). We ex- 
tracted the Chandra lightcurves and spectra using programs written in IDL 
version 5.4. The extraction areas for source and background spectrum are those 
defined in the Chandra User's Guide. 

The time of observation for the individual X-ray instruments is given in 
Table 1. We display the corresponding X-ray lightcurves in Fig. 1. The orbital 
phase has been computed from the ephemeris of Torres & Ribas (2001). First 
inspection reveals strong variability throughout the whole observation, including 
two large flares, and two 'high states' (i.e. extended phases of enhanced emis- 
sion) near the end of the observation. The secondary eclipse is clearly identified 
as a dip in the lightcurve close to orbital phase 0.5. Note, that the minimum of 
the X-ray lightcurve is not exactly centered on <E> = 0.5, but slightly offset to- 
wards earlier times. As we have observed simultaneously with two independent 
satellites a timing error is very unlikely. This shift may indicate an inhomoge- 
neous distribution of emitting material in the coronae of the YY Gem binary. 



Table 1 . Observing log for the XMM-Newton and Chandra observa- 
tions of the Castor system on Sep 29/30, 2000. 



Instrument 


UT 


JD - 2451817 


Expo 




Start 


Stop 


Start 


Stop 


[ksec] 


XMM-Newton 


EPIC-pn 


18:56 


08:57 


0.2888 


0.8730 


50.47 


EPIC-MOS 


18:15 


08:51 


0.2604 


0.8689 


52.57 


RGS 


18:07 


09:27 


0.2542 


0.8937 


55.26 


Chandra 


LETGS 


21:30 


13:54 


0.3958 


1.0792 


59.00 



The combination of XMM-Newton and Chandra allows to examine the X- 
ray spectrum of YYGem with intermediate (EPIC) and high (LETGS, RGS) 
resolution, and to compare the performance of the grating instruments on both 
satellites. The CCD spectra obtained with the EPIC are analysed in the XSPEC 
environment (version 11.0.1). 

3. The XMM-Newton EPIC spectrum 

We start with the analysis of the quiescent spectrum observed prior to the first 
large flare in Fig. 1 (JD 2451817.256 - JD 2451817.450). Following Giidel et al. 



Oh 

U 

3 



Oh 

o 



10 
8 
6 
4 
2 
1.0 
0.8 

0.6 



Orbital Phase 
0.0 0.2 0.4 0.6 



APS Conf. Ser. Style 
0.8 



1.0 





XMM EPIC-pn - 







i 0- 4 ^ 
* 0.2^ 
0.0 L_ 



0.2 0.4 0.6 0.8 

HJD- 2451817 



CHANDRA 



1.0 



Figure 1. X-ray lightcurve of YYGem observed by XMM-Newton 
(top panel) and Chandra (bottom panel) on Sep 29/30, 2000. The 
Chandra lightcurve has been extracted from the zeroth order image. 
For XMM-Newton we show the lightcurve observed by EPIC-pn. The 
orbital ephemeris is from Torres & Ribas (2001). 



(2001) we represent the quiescent EPIC spectrum of YYGem by a 3-temperature 
(3-T) model for thermal emission from an optically thin plasma (VMEKAL). In 
order to better constrain the spectral model we analyse the spectra from the 
pn and the two MOS detectors simultaneously. For the joint modeling of the 
spectrum from these three instruments we add a constant normalization factor to 
make up for uncertainties in the absolute calibration of the detectors. The EPIC 
spectrum for the pre-flare phase is shown in Fig. 2, and the best fit parameters 
from the 3-T model are summarized in Table 2. 

We use this spectrum as a baseline for time-resolved spectroscopy. The 
EPIC lightcurve of YYGem is split in a total of 15 time intervals (listed in 
Table 3) representing different activity levels of the star, and the spectrum of 
each phase is modeled by a 3-T model. As the integrated light from the quiescent 
corona should be visible at all times we hold all temperatures and abundances 
fixed on the values given in Table 2, and vary only the emission measure. In 
some of the time segments, namely for the post-eclipse feature and during the 
large flares, the 3-T model does not provide an adequate description of the 
EPIC spectrum: A high energy excess stands out in the residuals suggesting 
the presence of higher temperature material in addition to the emission from 
the quiescent corona. Adding a fourth VMEKAL component does not lead to a 
significant improvement. Only a 5-T model represents the data well (Xred ~ -0 
during the phases of most intense emission. For the modeling of these time 
intervals we have fixed spectral components # 1 - 3 on their quiescent values 
(see Table 2) . All abundances of components # 4 and # 5 have been held fixed on 
solar values because the statistics do not allow to constrain further parameters. 
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Table 2. Spectral parameters for the quiescent state of YYGem (de- 
rived from t < JD 2451817.450). Normalization constants for cross- 
calibration of the three instruments (pn, MOS 1, MOS2) are: N pn = 1 
(fixed), iV mosl = 1.0llg;°i, and N mos2 = 1.03l°;°J. 



fcTi 


feT 2 


kT 3 


[keV 




n 91 +0.05 
u - zi -0.07 


n R4+ - 01 

u - D4 -0.02 


i 7 n+0.30 








EM 2 


i?M 3 


[W bi cm" 3 




2 24 +L87 

z -^-0.60 




o OO+1.05 
z - oo -0.66 









Mg 


Si 






o.64±8:S 


n 27 +0 - 12 


u - 4 ' -0.07 






s 


Fe 


Ni 


Xrcd ( dof ) 


50 +u - 2V 

U.OU_ 20 


23 +u - 04 

u -^°-0.03 


n nn +u - 28 

u - uu -0.00 


1.24 (658) 



Table 3. Time intervals selected for a systematic investigation of the 
evolution of spectral parameters throughout the XMM-Newton EPIC 
observation from 29/30 Sep 2000. 



Start 


Stop 


Remarks 


Interval 


[JD - 


2451817] 






0.290 


0.450 


pre-flare quiescence 


h 


0.450 


0.475 


hump before flare 


ti 


0.475 


0.484 


rise flare 1 


H 


0.484 


0.493 


decay (a) flare 1 


U 


0.493 


0.510 


decay (b) flare 1 




0.510 


0.525 


mini-flare 


h 


0.525 


0.630 


post-flare quiescence 


h 


0.630 


0.662 


secondary eclipse (1st half) 


t 8 


0.662 


0.688 


secondary eclipse (2nd half) 




0.688 


0.710 


post-eclipse feature (a) 


ho 


0.710 


0.735 


post-eclipse feature (b) 


hi 


0.735 


0.775 


post-eclipse feature (c) 


tl2 


0.775 


0.790 


rise flare 2 


h 3 


0.790 


0.805 


decay flare 2 




0.805 


0.869 


'high state' 


tl5 
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Figure 2. XMM-Newton EPIC spectrum, left - Quiescent spectrum 
(ii) of pn, MOS 1, and MOS 2 with x 2 residuals; right - Three different 
activity states seen with pn: Quiescence (ti), flare rise (£3), and 'high 
state' (tis); note the Fe K-shell emission at ~ 6.7keV during £15. 



The last time interval (£15; the 'high state') is an exception: The signal at 
high energies is larger than for all other time segments, and broad Fe K-shell 
emission is clearly visible (see Fig. 2). We find an acceptable solution in this 
case for § = 0.47±g;^. 

3.1. Temperature - Emission Measure Diagrams 

The evolution of temperature and emission measure puts important constraints 
on the dynamics during flare decays. In a one-dimensional hydro-dynamic ap- 
proach to model stellar flares developed by Reale et al. (1993) the duration of 
the heating determines the slope in the lgT — lg (V EM)— diagram. We have de- 
rived lgT - lg ( v / EM)-dia grams for the spectral components of the 5-T model 
that represent the heated plasma during the two large flares, i.e. VMEKAL 
components # 4 and # 5. Fig. 3 shows the evolution of the two large flares, 
both starting with the rise phase (time interval £3 and £13, respectively). Under 
the assumption that the flare emission is concentrated in a single loop the slope 
C observed during the decay phase can be used to obtain an estimate for the 
loop half-length L. This method has been calibrated for several instruments 
including EPIC-pn (F. Reale, priv. comm.). We apply the equivalent of Eq. 2 
from Reale et al. (1997) to derive L from the slope £, the observed temperature 
(T max = 39 MK), and the decay constant of the lightcurve (t\ c = 16 ± lmin). 
The resulting loop length is L ~ 2 • 10 9 cm. 



4. High-resolution Spectra: XMM-Newton RGS and Chandra LETGS 

A comparison of the time-averaged first order X-ray spectra of YY Gem as ob- 
served with LETGS and RGS is given in Fig. 4. We only show the region 
between 10 — 26 A, which contains the strongest lines. Line identifications are 
given on top of the diagram. The spectrum is given in units of cts/s/bin. Since 
the RGS and the LETGS observations overlap for about 75 % in time, the rela- 
tive strength of the lines measured by both instruments should be similar, with 
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Figure 3. Temperature - emission measure - diagram for the two 
flares in the Sep 2000 observation of YY Gem. left - VMEKAL compo- 
nent # 4, right - VMEKAL component # 5. The numbers next to the 
data points indicate the respective time intervals from Table 3. 



some dependence of the line strength on the binsize, and the absolute numbers 
demonstrate directly the difference in sensitivity between RGS and LETGS. 

The Lya line of H-like O VIII is by far the strongest line in the spectrum 
with the highest photon flux, i.e. taking account of the effective area. Next 
to a number of iron L-shell transitions we identify the He-like triplets of four 
elements: Si XIII, NelX, OVII, and NVI. The OVII triplet is the strongest 
triplet and the only one which is clearly resolved and not blended with other 
lines. A detailed investigation of the properties of the coronal plasma making 
use of line ratios will be presented by Stelzer et al., in prep. 



5. Summary 

The X-ray lightcurve of YYGem shows that the object was subject to strong 
variability including two large outbursts during the time of observation. The pa- 
rameters of a 3-T model for the quiescent emission are compatible with results 
from the analysis of an earlier XMM-Newton observation of YY Gem presented 
by Giidel et al. (2001). Time-resolved modeling of the EPIC spectrum reveals 
the presence of a high temperature plasma (kT max = 3.4 keV) in flares. Ac- 
cording to a one-dimensional hydrodynamic model the flare emission arises in a 
semi-circular loop with ~ 2 • 10 9 cm length. This approach is certainly a sim- 
plification of the real situation which does involve a multi-temperature plasma 
and possibly complex loop systems. Nevertheless, the hydrodynamic approach 
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Figure 4. Comparison of the simultaneous XMM-Newton RGS 1 and 
RGS 2, and Chandra LETGS spectrum of YYGem in the range A = 
10 — 26 A. Straight horizontal lines in the RGS spectrum represent data 
gaps due to CCD failure or gaps inbetween individual chips. 
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is important: while simple quasi-static modeling tends to reproduce large loops 
the method applied here demonstrates that the coronal structures are likely to 
be much smaller than the radii of both stars in the YY Gem system. 

The simultaneous observation of YY Gem with Chandra and XMM-Newton 
demonstrates the different sensitivity of these instruments. Each of the two 
RGS provides roughly the same count rate as the LETGS first order spectrum. 
The LETGS is more sensitive at short wavelengths (see e.g. the region around 
the NelX triplet), while the sensitivity of RGS is slightly higher towards longer 
wavelengths (e.g. near the OVIII Lya line). 
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